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Abstract 

We study the effect of spatially correlated noise on coherence resonance (CR) in a Watts-Strogatz 
small-world network of Fitz Hugh-Nagumo neurons, where the noise correlation decays exponen- 
tially with distance between neurons. It is found that CR is considerably improved just by a 
small fraction of long-range connections for an intermediate coupling strength. For other coupling 
strengths, an abrupt change in CR occurs following the drastic fracture of the clustered structures 
in the network. Our study shows that spatially correlated noise plays a significant role in the 
phenomenon of CR through enforcing the clustering of the network. 
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The response of nonlinear systems to noise has attracted large attention. Especially 
stochastic resonance (SR) has been studied very extensively during the last two decades due 
to a number of applications in many fields, from physical to biological systems □ □ □ □ □ 
The main result of SR, which is somewhat counterintuitive, shows that noise at a proper 
strength optimizes the response of a nonlinear system to a subthreshold periodic signal. 
An optimal strength of noise can induce the most coherent motion in the system. SR-like 
behavior of the coherent motion can be induced purely by noise even in the absence of an 

n 

external periodic signal for a system at a saddle- node bifurcation point [6] . This phenomenon 
has been called coherence resonance (CR) or autonomous SR. In general, various excitable 
systems such as the Fitz Hugh-Nagumo model |7(, the Plant model, the Hindermarsh-Rose 
model ^] and the Ho dgkin- Huxley model exhibit such noise induced coherent motion. 

Recently, SR and CR in coupled or extended systems have become an interesting is- 
sue, and some new features have been demonstrated, namely the noise-enhanced phase 
synchronization , the noise- induced spatiotemporal pattern formation ^j] and the noise- 



enhanced wave propagation 



Also, the phenomena called array enhanced stochastic 



have 



resonance (AESR) and array enhanced coherence resonance (AECR) 
drawn interests among researchers in recent years. It is now understood that in spatially 
extended systems (i) the topology of connecting structure and (ii) noise correlation among 
the elements are the significant ingredients on the collective behavior of the systems. Actu- 
ally, the connection topology of a variety of extended systems can be described by complex 
networks Especially many biological neural networks present clear clustered structure 
and sparsely long-range random connectivity [20]. On the other hand, spatially correlated 
noise has been considered relevant for biological systems. However, most previous studies 



have not dealt with these factors together |R 



13, 



18, 



2l|. In this study, we plan to add 



the element of spatially correlated noise in the system of excitable cells and to study the 
corresponding dynamics of CR. 

As a model, we consider a system of coupled excitable cells on "small-world" network, 
introduced by Watts and Strogatz [22| in the presence of a spatially correlated noise. Each 
cell is a Fitz Hugh-Nagumo (FHN) neuron which is a simple but representative model of 
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excitable neuron 23j] . N neurons in a ring lattice are diffusively coupled as following: 

%i Hi / j 9ij jp^j (1) 

3 

y\ = Xi + a + ^i, (2) 

where X{ is the fast voltage variable and yi is the slow recovery variable of ith neuron, e and 
a are a time scale and a bifurcation parameter respectively. We fix e = 0.01 and a = 1.03 for 
all N = 101 neurons. For \a\ > 1, a single FHN neuron has only a stable fixed point. While 
for \a\ < 1, a limit cycle occurs, gij is a coupling strength between two neurons % and j. 
If connected, they have the coupling strength = g, otherwise g^ = 0. The connectivity 
pattern can vary with parameter p, which measures the network randomness. & is spatially 
correlated noise with intensity D. 

The spatially correlated noise £ is generated by summing N Gaussian white noises ( with 
correlation function C . 

& = /sr r2 $Z ^+ kCk ( 3 ) 
V 2^keA u k keA 

A = {-4A,... ,-2,-1,0,1,2,. ..,4A}, (4) 

where Q is a Gaussian white noise with zero mean and correlation given by {Ci{t)Cj{t')) — 
D5ij5(t — £'); D denoting the noise intensity. And the correlation function is defined by 
Ck = exp(— 2k 2 /X 2 ). According to the above method, we can get the spatially correlated 
noise which obeys the correlation function with decay constant A as following: 

mMt')) = De W (-^f^) S(t - If), (5) 

where i and j denote the spatial positions of neurons in the ring lattice. Consequently \i —j\ 
represents a distance along the ring lattice not the connection topology. The systems are 
numerically integrated by the method of Fox et al. J2J] with the time step At = 0.002 t.u. 
(time units). 

As a quantitative observable of a neuron showing CR, a temporal coherence is calculated 

by 

S, = ^= x = l,2,...,JV, (6) 

where r« is the ensemble of time interval of inter spike and Si denotes temporal coherence 
factor of ith neuron. Here (-) t denote average over time. The coherence factor S of the 



system is computed by averaging Si over all N neurons. A larger S implies that the inter 
spike intervals of neurons are more uniform. 

A general feature of SR and CR is that there exists an optimal noise intensity at which 
the coherence factor is maximized. In the cases of AESR and AECR there exists an optimal 
coupling « additionally UUU- Our —1 stations verify the above 
results (see Fig. [TJ), that is, both the optimal noise intensity and the optimal coupling 
strength exist, moreover, regardless to the network randomness p and the correlation length 
A of the noise. For a very weak coupling g rs 10~ 2 ' 5 , each element is effectively independent. 
Therefore a coherence factor curve has the similar appearance to that of a single uncoupled 
neuron with the same parameters. As coupling strength increases, the impact of noise- 
induced firing events propagate through the connection topology. This cooperation of noise- 
induced individual and coupling-induced mutual excitation enhances a coherent motion in 
the coupled system. However, a strong coupling rather disturbs excitation of elements 
because those elements that are excited by noise are strongly attracted to a resting state 
by resting neighbors. For this reason, a stronger noise is needed to overcoming stronger 
coupling and to excite the neurons. In this case a global synchronization emerges due to 
strong coupling while the temporal coherence of the system is somewhat reduced due to 
strong noise. As seen in Fig. ^ the resonance curve shifts to right side gradually and its 
peak rises up first and then drops down according to the coupling strength. 

How do the connection topology and the spatial correlation of noise influence the coupled 
excitable neurons? To focus on the question the effect of structural changes of connectivity 
on the spatial correlation of noise must be explained in advance. For a regular network 
(p = 0), if A is large enough, each neuron interacts only with the ones exposed to the 
correlated noise. As p increases each neuron is able to interact with a distant one which 
is exposed to an uncorrelated noise. Consequently for a proper length of A, increasing of 
p reduces the correlation of the noise and simultaneously increases the small-world effect 
between coupled neurons. However, when A = every neuron is already exposed to a 
totally uncorrelated noise, i.e. local noise, so varying p does not alter the correlation of 
noise. In this case the variation of p influences the small-world effect only. 

Now let us figure out the peculiar result according to p and A. In Fig. ^the maximum 
value of each curve is called a maximal coherence factor S m and the noise intensity at S m 
is called an optimal noise D opt . If the coupling strength is not very strong, S m takes a 
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larger value for completely random network (p = 1) than for regular network (p = 0). On 
the other hand, as the coupling strength increases further S m shows a rapid decline for the 
more random network. For large p the interaction among neurons becomes more fast due 
to long-range connections which amplify the effect of coupling for the entire range of g. 
Therefore, when the coupling is not very strong the long-range connections enhance further 
the coherent motion of the system. However, for the very strong coupling strength, the 
stronger D opt is required due to the reinforced synchronization by long-range connections. 
Hence S m decreases rapidly for large p. The result of Fig. IHa)-(b) corresponds to the 
study of Ref. Q|. Figure ^c)-(d) show the results for the cases of non-zero A. The overall 
tendency does not change except that the coherent motion is depressed when compared to 
the case with an uncorrelated noise. Generally the correlation of noise makes the firing 
events of neurons correlated so the effect of mutual excitation is diminished. As a result 
a coherent motion is depressed r], 21]. Interestingly, we found that when p = and 



A = 2, D opt has a relatively low value even for a large coupling strength. When a neuron is 
excited by a noise, its neighboring neurons are excited by the coherent noise most probably. 
Consequently, spatially correlated noise enables neurons to generate firing events in the 
smaller noise intensity even for a relatively strong coupling. However, because the correlation 
of noise disappears mostly again for p near 1, D opt s return to those values of the case for 
p = 1 and A = 0. 

We then systematically examined the maximal coherence factor S m as a function of 
network randomness p for different values of A (see Fig. EJ). We fix coupling strengths at 
five values; very weak (g pa 1CT 2 ' 5 ), rather weak (g pa 1CT 2 ' 25 ), optimal (g ~ 10~ ), rather 
strong (g pa 10 25 ), and very strong (g pa 10~ ' 75 ). For some cases at A = (see Fig. 
I^c)-(d)), S m rises rapidly as p increases when the coupling strength is not very strong. And 
S m rather decreases as p increases for a stronger cou plin g g pa 10~ 75 (see Fig. Efe)). The 
similar behaviors have also been observed in the SR 17J. 

Some interesting features are also found in a certain region of parameter space. Contrary 



to an expectation from the existing studies |lfl|21|], the spatial correlation of noise enhances 
a coherent motion. If the coupling strength is very weak and also if A is zero, each element 
behaves independently. However, for A > 0, S m curve is elevated slightly according to A 
(see Fig. |2Ja)) since the spatially correlated noise assists partially synchronized excitation 
among the neighboring neurons even for a very weak coupling. As p increases, the effect of 
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the correlated noise disappears and the partial synchronization vanishes too. Therefore, for 
each value of A, S m decreases with increasing p and finally reaches the value at A = 0. For 
other values of coupling strength, however, coherence deteriorates in general by the spatial 
correlation of noise. 

All other S m curves in Fig. |2] are divided roughly into two categories depending on 
whether S m curve monotonically increases or not. In the case of rather weak coupling the 
effect of long-range coupling is clearly observed. When A is large for this case, we can define 
the transition point of p around p c = 0.1 above which S m grows significantly. In Fig. El the 
clustering coefficient C remains practically unchanged (clustered structure mostly remains) 
for p < p c while the characteristic path length L drops sufficiently As p passes over p c 
C drops rapidly while L rarely changes. It indicates that introducing a few long-range 
connections is enough to decrease L sufficiently and additional long-range connections for 
p > p c affect only to fracture the clustered structure. Because of the clustered structure 
of neurons consolidated by a spatially correlated noise, a few long-range connections do 
not affect the elevation of the S m at large A for a rather weak coupling. As p is increased 
beyond p c , the cluster begins to be fractured and then the effect of noise correlation vanishes 
considerably. This makes S m , which is definitely dropped due to the correlated noise at 
small p, grow rapidly to the value it has when A = at large p. Since the coupling term 
in Eq. ^ affects the dynamics for a larger coupling strength, the change of connectivity 
becomes more crucial. By every rewiring to the long-range connections the coherent motion 
is steadily enhanced as p increases for entire A. This ascending behavior of S m due to the a 
few long-range connection is maintained to a rather strong coupling (see Fig. |2fc) and (d)). 

When we consider both the noise correlation and the connection topology, we observe 
another peculiar behavior of S m . S m has a maximum value at an intermediate p around 
p c for a rather strong coupling when A > (see Fig. Efd)). In the previous studies S m 
either increases or decreases monotonically For this to make sense, we need to 

follow the behavior of D opt for each situation. When A = the optimal noise intensity D opt 
corresponding to S m rarely changes for the entire range of p (see Fig. Ufa)). In case of A > 0, 
D opt remains constant at lower level than in case of A = until p « p c and after that point 
gradually chases the value of A = (totally uncorrelated noise). As has been mentioned 
above, the correlated noise enables resonance curve to occur in smaller noise level especially 
for relatively strong coupling. This indicates that the effect of correlation of noise between 
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inter-neurons drastically vanishes as clustered structure is fractured. As p increases further 
from p c for A > 0, resonance curve shifts to higher noises. This is why the value of S m does 
not increase and even decreases after p > p c . 

Finally for the very strong coupling (see Fig. I2fe)), S m changes a little until p = p c 
then decreases substantially for p > p c . This depression is more distinct as A increases. 
A few long-range connections do not affect S m much, because of the clustered structure 
consolidated by strong coupling and moreover correlated noise especially for large A. The 
behavior of D opt as a function of p for various A is very similar to the case of rather strong 
coupling as seen in Fig. Hfb). S m should decrease much as A increases when p is near zero. 
However the decrease is not big owing to the effect of the reduced D opt . This tendency also 
appears in rather strong coupling case. As a whole, for p > p c , the movement of D opt toward 
larger value and a number of long-range connections with very strong coupling drastically 
reduce the temporal coherence of the system. 

In summary, we have investigated the effect of spatially correlated noise (correlation 
length A) in the presence of various connection topology (network randomness parameter 
< p < 1) of FHN neural network. This study reproduces most of the general feature of 
AESR and AECR in the entire range of A and p. In addition, we could get some novel 
features for A > 0. When the coupling is rather strong, an optimum value p m p c clearly 
emerges where a maximal coherence resonance appears. For p increasing from to p c , the 
maximal coherence factor rarely changes but, beyond p ~ p c , it grows dramatically either 
up for a weak coupling or down for a strong coupling. For p beyond p c , it is observed that 
the clustered structures of the neurons are mostly fractured out due to many long-range 
connections and, as a result, the noise correlation of inter-neurons is diminished quickly. It 
is believed to be the reason for such an abrupt change of S m near p w p c . These results show 
that the spatially correlated noise enforces the role of the clustered structure to the system. 
Therefore, the effect of a few long-range connections is ignored by the enhanced clusters of 
the neurons for a weak as well as for a strong coupling. Nevertheless, for an optimal and 
for a rather strong coupling, coherence resonance is still considerably enhanced by a small 
portion of long-range connections. 
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FIG. 1: The coherence factor S versus noise intensity D for several values of coupling strength on 
(a) regular network (p = 0) and (b) completely random network (p = 1) with spatial correlation 
length of noise A = 0. (c), (d) are the same curves for p = and p = 1 with A = 2 respectively. 
Average number of neighbors k = 6. Number of elements N = 101. Coupling strength g varies 



from g = 10 to g = 10 with 0.25 step of exponent. Base curves represent for g = 10 



-2.5 



As coupling strength increases, the peak of S shifts to a stronger noise. 
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FIG. 2: The maximal coherence factor S m as a function of network randomness p for different 
values of A, when coupling strength is (a) very weak g = 1CT 2 - 5 , (b) rather weak g = 1CT 2 - 25 , (c) 
optimal g = 1CT L75 , (d) rather strong g = 1CT L25 , and (e) very strong g = 1CT - 75 . A is varied 
from A = (local noise) to A = 6 with a step 2. S m curve goes down with increasing A except in 
(a). 
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FIG. 3: Clustering coefficient C (open circle) and characteristic path length L (solid circle) as 
a function of p for Watts- Strogatz small-world networks with N = 101 and k = 6. They are 
normalized by each value at p = 0. 
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FIG. 4: Optimal noise intensity D pot as a function of p for various A including (diamond), 2 
(circle), 4 (triangle) and 6 (rectangle), (a) g = 1(T L25 and (b) g = 10~ a75 . 
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